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Estimation of parameters is a pivotal task throughout science and technology. In quantum mechanics, quan- 
tum Cramer-Rao bound provides a fundamental limit of precision allowed to achieve under quantum theory. For 
closed quantum systems, it has been shown how the estimation precision depends on the underlying dynam- 
ics. Here, we propose a general, alternative formulation for an estimation scenario in open quantum systems, 
aiming to relate the precision more directly to properties of the underlying dynamics. Specifically, we derive a 
Cramer-Rao bound for a fairly large class of open system dynamics, which is governed by a Lindbladian master 
equation. We illustrate the utility of this scenario through an example. 
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Introduction. — Metrology and parameter estimation lie in 
the heart of science, and are prevalent in any aspect of tech- 
nology. The basic task of identification or estimation of a set 
of unknown parameters essentially requires an inference from 
a pool of observed data about the parameters or the system 
to which they are attributed. As errors and imperfections are 
unavoidable in practice, increasing accuracy of the underlying 
tasks of data acquisition and inference — hence improving the 
quality of estimation — is an important goal of metrology (TJ. 
Improving quality of measurement instruments and removing 
potential sources of systematic errors aside, statistics provides 
useful suggestions for enhancing metrology, such as increas- 
ing data size and repeated measurements on an ensemble of 
N 'probe' systems. Additionally (and more interestingly), 
the underlying physics of the system of interest may also dic- 
tate some restrictions or bounds on the ultimate achievable ac- 
curacy (usually described through a 'Cramer-Rao inequality' 
[2 1), or even may offer new possibilities to exploit. 

In quantum mechanics, measurements act differently than 
in classical systems. In addition, interactions with an environ- 
ment or other systems as well as (quantum) correlations can 
each affect observed data |3;|, hence introduce new playing 
factors in estimation theory. For example, it has been shown 
that entanglement in a probe ensemble can be exploited to the 
advantage of a quantum estimation task [4], so that it enables 
the estimation error of 0(1 /N) (the "Heisenberg limit"), in 
contrast to the classical statistical limit of 0(1/ yN) (the 
"shot-noise limit"). Alternatively, enabling fc-body (k ^ 2) 
interactions among quantum probe systems has been shown 
to allow an error of 0(1/ V N k ) 0; or, it has been argued 
that application of a suitable entangling operator may even 
offer an error as small as 0(2~ N ) J6| (beyond the Heisen- 
berg limit). Moreover, nonclassicality has been examined as 
a potential resource for increasing the metrology resolution in 
quantum optics [7| (for a general framework of resource anal- 
ysis, see, e.g., Ref. El). It thus seems natural to expect that 
some properties of quantum systems can be employed as a 
useful "resource" for quantum metrology. 

In open quantum systems, due to interaction with an en- 
vironment, the underlying dynamics becomes 'noisy.' As 
a result, formulation and analysis of quantum estimation 
also becomes more involved (9] [TO]. In general, dy- 



namics of an open system can be described as Qs(i) — 
Tr E [UsE(t, t )gsE(to)ul E (t, t )}, where q S e is the state of 
the systems and environment (SE), and UsE(t, to) is the cor- 
responding unitary evolution ifTTl [T2l . Thereby one can ar- 
gue that in general there may exist a flow of information be- 
tween the system and the environment |[T3ll . Under some con- 
ditions, this dynamics can feature quantum Markovian or non- 
Markovian properties [ 14 1. The former case typically appears 
when the environment has a small decoherence time during 
which correlations typically disappear, whereas in the latter 
correlations (both classical and/or quantum lfl6l ) with the en- 
vironment would form and persist. Such correlations are in 
practice inevitable, which necessitates investigation of 'noisy' 
quantum metrology J5] [17] [THIl . Additionally, such correla- 
tions may also offer new resources for enhancing estimation 
tasks. A general analysis of the latter possibility is still lacking 
(see Ref. |[T9l for an example along this line). 

In this article, we first lay out a fairly general formalism 
for open quantum system metrology. This (re)formulation of 
the problem (cf. Ref. [9]) has this advantage that here pre- 
cision of estimation is more directly related to the underly- 
ing dynamics; besides, it is in some sense analogous to the 
closed system formulation. In particular, we derive a quan- 
tum Cramer-Rao bound (QCRB) for open system dynamics 
generated through dynamical map with semigroup property. 
We next illustrate this setting through an example. This ex- 
ample shows that how induced correlations of probe quantum 
systems through a common environment may offer relatively 
higher precision for an estimation task in a sense akin to what 
manybody interactions enable. 

Open system dynamics. — Under some specific conditions, 
the dynamical equation describing the state of an open sys- 
tem gg [defined on a Hilbert space Hs, i.e., Qs € S(Hs), 
where S(Hs) is the space of all linear operators acting 
on H] reduces to d T gs(r) — C t [qs(t)], or equivalently 

Qs(t) = Te^ T o £t,cIt [f?s( T o)]> in which T denotes time- 
ordering, and £ T [o] = -i[H s (T),o] + ]T fc T}k(r)(A k (r) o 
A\(t) - (l/2){A\(T)A k (T),o}) (for some set of operators 
{Ak(r)}) is the (Lindbladian) generator of the dynamical 
map, with Hs (r) being the system Hamiltonian up to a Lamb 
shift term (we omit subscript S henceforth). We have also 
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assumed H = 1. In (time-dependent) Markovian evolutions, 
we have J7fc(r) ^ Vfc, r; while if some rjk becomes neg- 
ative for some intervals, the associated dynamics would be 
non-Markovian ifTTI [TH [141 [151 . 

Let us assume that a set of unknown parameters x = 
(xi, . . . , xe) are to be estimated in a quantum system subject 
to interaction with an environment. For simplicity, here we 
only consider the single parameter case (x). In the closed- 
system scenario, this parameter is usually assumed to enter 
into the dynamics as a linear coupling in the Hamiltonian 
H(x) = xH acting on some known initial state. In the open- 
system scenario, similarly the devised dynamics would in gen- 
eral depend on x, whence the dynamics of the system is de- 
scribed by d T g(x,t) — C T (x)[g(x, t)]. For our analysis in 
this article, as we argue later, vectorization is an appropri- 
ate language. Vectorization of this equation yields d T | g)) = 
C T (x)\g)), where C T is the matrix representation of C T ll20l . 
Defining the normalized pure state g = \g)) ((o|/Tr[p 2 ] £ 



have 



and assuming that C T does not depend on time, we 



q(x,t) 



,r£t (x) 



/Tr[e r£(x) p(0)e r£t(x) ] 



(1) 



The initial preparation g(0) may itself depend on x, but here 
we do not assume such generality. 

QCRB for open system estimation. — Given a data set D = 
{7i} constituted from some measurement outcomes ji over 
N (identical) probe systems, an estimator x est (D) is chosen 
for the true value x. By repeating this scenario M times and 
averaging, the precision of the estimation of x, evaluated by 
5x = ^/Var(x), is then fundamentally limited by the QCRB 



5x ^ l/VMyV(Q)(x; N). 



(2) 



Here, Var(x) is the variance of any unbiased estimator x est (D) 
(for which, by definition, (x est ) = x, with (o) denoting the av- 
erage with respect to the underlying quantum probability dis- 
tribution), and T^ix-jN) is the so-called "quantum Fisher 
information" (QFI) ifTTl I2T1 l22l . By assuming the state of 
each N -probe set to be g( N \x,r) (hereafter we omit super- 
script N for brevity) and defining the corresponding symmet- 
ric logarithmic derivative L g through <9 x p = {L e g + gL B )/2, 
the QFI is defined as T^(x, r; N) = Tr[g(x, t)L 2 ,1 

We remind that in closed systems, noting g(x, r) = 
e -t-r-ff(x) £(o)e lTj? ( x ), the spectral decomposition g — 
< £r i \r i )(r i \, and L e = 2Y, ij (ri\d x Q\r j )/(r i + r j )\r i )(r j \ 
provides a relatively straightforward relation between J 7 ^) 
and the interaction H. In particular, when H(x) = xH and g 
is pure, one obtains 



jr(Q) =4t 2 Cov (H,H) 



(3) 



(with equality replaced with ^ for mixed g), where 
Cov g (X, Y) = (XY) g - (X) g (Y) e is the covariance func- 
tion of a pair of observables X and Y, which here is the very 



resulting relation 

5x > l/(2r\/My / Cov e (H,H)) = l/(2rV / MAH), (4) 

where AH = \/A 2 H, is more in the spirit of an uncertainty- 
like relation BTTl . and shows explicitly how the precision is 
dictated by the interaction. In open-system cases, however, 
deriving similar, straightforward relations is hardly possible 
(except for some special cases lfT3l ) since, e.g., calculating 
L e(x,r) = 2/ °°e~ se ( x ' r )d x e(x,T)e- se ( x ' r )ds is involved. 
Thus it is difficult to capture how interaction with an envi- 
ronment affects the QFI and the precision. To partially alle- 
viate this issue, here we follow an alternative approach work- 
ing with the vectorized state g instead, which enables a bound 
somewhat akin to Eq. Q — with H replaced with C. Although 
an obvious expense to be paid in going from g to g is the (ar- 
tificial) nonlinearity of relations in terms of p, we show that 
this formalism still retains its relative utility in providing use- 
ful bounds. Let us briefly remark that the mentioned nonlin- 
earity problem could be lifted if we instead worked with the 
purification w of g = ww\ but this does not really offer an 
advantage over the original picture in leading to a straightfor- 
ward uncertainly-like relation. 

Now from the symmetric logarithmic derivative Lg — 
2d x g, one can define an associated QFI ^F^> by replacing 
(g,L g ) — > (g,Lg) in J 7 ^. After some straightforward al- 
gebra, using the dynamical equation Eq. ([TJ, and assuming a 
linear x-dependence as C(x) = xL, it can be seen that 



7<Q) =4r 2 CoVg(L t ,L). 



(5) 



This relation is analogous to Eq. Q, where instead of the 
Hamiltonian we have the Lindbladian generator of the open 
dynamics. 

The QFI has a natural interpretation. Recall that 
indeed emerges from the optimization of the Fisher informa- 
tion over all possible quantum measurements on the system 
g £ S{%) ||2TI . Similarly then, J-^ is obtained if any quan- 
tum measurement on the 'system' g £ S(H® 2 ) is allowed. 
Note, however, that a natural extension of the measurements 
in H to H® 2 does not necessarily translate into most general 
measurements there. For example, a complete set of measure- 
ment {n^} (with the properties n, ~>t and J^IL; = 1-^), 
when extended simply as = |IL;)) do not constitute a 
complete set in the sense that £V II; 7^ 1%®2 in general. 

Let us see how compares with . First we re- 

mark that, from vectorizing the very definition of the symmet- 
ric logarithmic derivative, we have Lg — L e <E> 1 + 1 
d x In Tr[p 2 ]. This in turn yields the following expression: 



■ (Tr[gL e gL e ] + Tr[p 2 L 2 ] - 



quantum standard deviation A H (with 



Tr[po]). The 



(6) 

This form is not directly related to . However, using the 
relation Tr[XF] > A min (X )Tr[Y] (valid for any pair of pos- 
itive matrices X and Y\ here A m i n (X) denotes the minimum 
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eigenvalue of X) and Tr[g 2 ] ^ 1, we obtain 

•F (Q) > 4A min (e)J-K» - F(f>), 



Probes 



Common Bath 



(7) 



where F(g) represents the last term in Eq. d6|l. This bound 
would be trivial where during the evolution the state of the 
system becomes pure. An example of such cases is when the 
dynamics is unitary and the initial state is pure (|W(x, r)) = 
U(x, t)|^(0)), with unitary U). For such cases a signifi- 
cant simplification occurs due to (^(x, r)|L*|W(x, r)) = 0, 
whence Eq. ^ reduces to = 2J 7 ^\ Therefore, in gen- 
eral we have 



1/J- (Q) > kj{T 



(Q) 



F), 



(8) 



where for nonunitary evolutions in which the state of the sys- 
tem never becomes pure we have k = 4A m i n (p), and for uni- 
tary evolutions with a pure initial state we have k — 2, F = 0, 
and the inequality is replaced with equality. In addition, for 
unitary evolutions with a mixed initial state g(0), we have 



4r 2 



Tr[p 2 (0)] 



Tr[ e 2 (x,T)H 2 -(g(x,T)H) 2 ]. (9) 



This result provides another framework for open-system es- 
timation, alternative to the approach of Ref. [9|. A natural 
advantage is that here the QFI is more directly related to the 
generator of the dynamics, and is straightforward to compute. 

Example. — Let us assume that there are N probe particles 
each of which only interacts with a common bath such that the 
interactions induce all possible fc-body terms (Fig. [T| l23l in 
the Lindbladian as follow: 



L[o] = £ a n 



C N , k o, (10) 



where er^ are all the same Pauli matrix (e.g., a z ), subscript i 
is the particle index, and the factor CV.fc = (^) counts the 
number of A; -body operators. 

We choose the initial state of the whole A -probe system to 
be the maximally entangled pure state g(0) = |4 r )(^ r | (i.e., 
|f>(0)» = where |*) = (\\ M )® N - \\ m )® N )/y/2, 

and A m (Am) is the smallest (largest) eigenvalue of a. For odd 
k, a h a i2 ■■■<Ti k ®<T il a i2 ---a ih (\y)®\y*)) = I*- 1 ) I*- 1 *), 
where I*- 1 ) = {\\ M )® N + \\ m )® N )/V2- It is straightfor- 
ward to see that 

L|*)|**) - C N , k (\* x )\* x *) - |*)|**)), (11) 

whence (|^ ) l^ 1- ) — |^)|^*})/v2 becomes an eigenvector 
of L corresponding to the eigenvalue — 2Cjv fc. Hence 



^)(x,t) 



16t 2 C 2 fc( 



-4C N ,fcTx 



[g-4Cjv,fcrx _ ^ e -2C]v,fcTx _|_ gji 



(12) 



An immediate implication of this relation is that for small val- 
ues of the x parameter, noting Cn ,k = 0(N k ) for large As, a 
polynomial precision in the estimation can be achieved. 



"Pi 


H BPi 




"p 2 


H BP 2 




■ 








H B P N 











FIG. 1. N probes, initially well isolated from each other, all interact 
with a common bath through two-body interactions HsPi- Here Hp i 
and Hb are the free Hamiltonians of probe i and the bath, respec- 
tively. These two-body interactions may induce a manybody quan- 
tum correlation among the probes |23 1. 



Summary. — Here we have outlined a fairly general for- 
malism for open quantum system metrology. In this 
(re)formulation of open quantum system metrology precision 
of estimation is more directly related to the underlying dy- 
namics, in some sense similar to the closed system formula- 
tion. The core of this formalism is the derivation of a quan- 
tum Cramer-Rao bound for open system dynamics generated 
through dynamical map with the semigroup property. This 
setting then was illustrated with an example. This example 
has implied that it may be possible to exploit induced correla- 
tions of probe quantum systems through a common environ- 
ment in order to achieve a relatively higher precision for an 
estimation task. 

Our formalism may introduce methods for utilizing some 
of the resources offered in open quantum dynamics, such as 
induced manybody correlations and memory, to hopefully en- 
hance a quantum estimation task in the presence of 'noise'. 
This in turn can open up possibilities for applications in , e.g., 
quantum sensing ifTTl |24H and quantum control of optome- 
chanical devices for advanced technologies [25|. 
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